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Abstract 

We discuss the contribution of ultraviolet (UV) renormalons in QCD to two-point func¬ 
tions of quark current operators. This explicitly includes effects due to the exchange of 
one renormalon chain as well as two chains. It is shown that, when the external euclidean 
momentum of the two-point functions becomes smaller than the scale associated with 
the Landau singularity of the QCD one-loop running coupling constant, the positions of 
the UV renormalons in the Borel plane become true singularities in the integration range 
of the Borel transform. This introduces ambiguities in the evaluation of the corresponding 
two-point functions. The ambiguities associated with the leading UV renormalon singu¬ 
larity are of the same type as the contribution due to the inclusion of dimension d = 6 
local operators in a low-energy effective Lagrangian valid at scales smaller than A^. We 
then discuss the inclusion of an infinite number of renormalon chains and argue that the 
previous ambiguity hints at a plausible approximation scheme for low-energy QCD, re¬ 
sulting in an effective Lagrangian similar to the one of the extended Nambu-Jona-Lasinio 
(ENJL) model of QCD at large N c . 


1 Work partially supported by research project CICYT-AEN95-0882. 



1 Introduction. 


The origin of renormalons in quantum field theory goes back to early work by ’t Hooft Q, 
Lautrup Q and Parisi Q . They made the observation that in renormalizable theories like for 
example quantum electrodynamics (QED), there exists a class of Feynman diagrams which 
give rise to a characteristic pattern in the coefficients of the large-order terms in the per¬ 
turbation series; these coefficients have an n!-growth with the same sign indicating that 
the corresponding series are not Borel summable. [] In QED this growth originates in the 
large momentum integration region of virtual photons dressed with vacuum polarization in¬ 
sertions which leads to singularities in the associated Borel plane; the so-called UV renor¬ 
malons j|,^,[5| |7|]. In quantum chromodynamics (QCD), it is the low momentum integration 
region of virtual gluons dressed with running couplings what is the source of the non integrable 
singularities in the Borel plane; the so-called IR renormalons jl],|8],[9|] . The appearance of these 
singularities is perhaps not surprising since in the extreme kinematic regimes in question it 
is not expected that these theories are well described by a simple perturbative expansion in 
the coupling constant. 

There are two good reasons to concentrate on the class of Feynman graphs which are at 
the origin of the renormalon singularities. One is the fact that in massless theories, like QCD 
with massless quarks, the breaking of conformal invariance is encoded in the /3-function. If 
one wants to find hints from perturbation theory on the possible origin of scales it seems ap¬ 
propriate to focus one’s attention on the summation of the infinite subset of graphs associated 
with the renormalization of the coupling constant. The other reason is the opportunity that 
is offered to explore issues which have to do with the analyticity (or rather, lack of it) in the 
coupling constant. 

The study of renormalon properties in gauge theories is at present an active field of re¬ 
search. The role of IR renormalons in the operator product expansion of two-point functions 
and their relationship with non-perturbative inverse power corrections has been extensively 
discussed in the literature 10- This and further discussions which originate in the work 
of ref. jn]], has led to a new point of view [|l^] concerning renormalons in QCD. The focus 
is now on the possibility that their systematic study in a given hadronic process might sug¬ 
gest generic non-perturbative effects of a universal nature. The basic idea is that genuinely 
non-perturbative effects ought to cure [Aany disease which appears when perturbation the¬ 
ory results are analytically continued |l9|]. This applies to the case of IR renormalons and, 


as recently suggested by Vainshtein and Zakharov(V-Z) ||d], [21j, perhaps to the much less 
explored role of UV renormalons in QCD as well. 


In QED (or for that matter in any renormalizable theory which is not asymptotically 
free) the UV renormalons cause a real obstruction to defining the theory in the ultraviolet 
from a resununation of the perturbative series. In earlier work by Parisi on </> 4 theories || 
it was argued that it ought to be possible to mimic the contribution from large momenta 
(namely UV renormalons) by means of insertions of local dimension six and higher opera¬ 
tors @ . Specifically, local here means that all the physics at momentum larger than a certain 
scale k is encoded in higher dimensional composite operators that are local on a scale of 
l/k. In other words, these composite operators are products of fields evaluated at the same 
point but suppressed by inverse powers of the cut-off momentum k. On these grounds, one 

2 For a comprehensive review of the subject and a collection of articles previous to 1990 see ref. 0. 
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would expect that by adding carefully-adjusted higher-dimensional operators to the initial 
Lagrangian one ought to be able to remove all the UV renormalons from QED to all orders 
in the coupling constant j22|. The intuitive reason is that QED can always be thought of 
as the low-energy limit of a larger asymptotically free theory (like e.g. SU( 5)) for which, 
in principle, UV renormalons are not a problem. Of course this means adding an infinite 
tower of higher dimensional operators which will make QED an effective nonrenormalizable 
theory, in agreement with its alleged triviality |23|]. Notice that the previous argument is 
entirely perturbative. In other words, the higher dimensional operators in question are still 
considered as irrelevant. If, on the other hand, a few of them turned out to be relevant, it 
might be possible to truncate the list of operators to just the relevant ones. This is actually 


what happens in the quenched ladder approximation of QED [ 241 where four-fermion opera¬ 
tors of the Nambu-Jona-Lasinio type [^5| (NJL) turn out to become relevant operators and 
consequently are kept in the Lagrangian on an equal footing as the four dimensional ones. 

In this connection we wish to point out that there has also been an ongoing struggle in the 
lattice community in trying to clarify if an abelian theory in the strong coupling regime has 
anything to do with the dynamics of models a la Nambu-Jona-Lasinio. Although this point 
is still unclear |26|] , there seems to be some evidence that four-fermion operators may play an 
indispensable role for understanding the QED non-perturbative (i.e. ultraviolet) dynamics. 


In QCD, and provided that the external momenta of a given Green’s function are not 
in regions of exceptional momenta [[l?]], the UV renormalons are not an obstruction to inte¬ 
grate virtual Euclidean momenta. Therefore one does not expect them to be at the origin 
of fundamental ambiguities. In their recent work |2(], 211, V-Z have found however that the 
contribution to the Adler function from the leading UV renormalon coming from the ex¬ 
change of two chains of vacuum polarization self-energy loops in an Abelian-like model was, 
contrary to naive expectations, dominant over the contribution coming from a single chain. 
Furthermore, they also argued that the contribution from more and more chains should be 
equally important and that, consequently, the actual value of the residue of the leading UV 
renormalon was ill defined. They then concluded that, at the phenomenological level, this 
could be taken as an indication that a new type of ambiguity may appear , which in the 
particular case of a two-point function with Euclidean momentum Q 2 , shows up as possible 
l/Q 2 -like contributions. 

Interestingly enough, the leading UV behaviour found by V-Z was shown to originate in 
the contribution from the insertion of d = 6 four-fermion operators. The appearance of these 


four-fermion operators, much like those of the Nambu-Jona-Lasinio model [25], is rather 
intriguing. It is known that there are extensions of this model which, when taken as models 
of QCD in the large-V c limit at intermediate scales < CJ(lGeV), are rather successful |28|, 29| 
in predicting low-energy physics (like for instance the Li coupling constants of the Lagrangian 
of Gasser and Leutwyler [30|.) However, the possible connection between these models and 
QCD has remained so far a mystery. 


There are some alternative routes to NJL-type models which have also been suggested to 
describe low-energy QCD |n],|32|]. Although not yet comparable at the level of phenomeno¬ 
logical success, they are nevertheless interesting in the sense that they represent relatively 
small departures from perturbative QCD and, hence, they ensure that at least in some limit 
they are likely to be related to QCD. Needless to say, the tough problem they face is to 
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show that the departures from perturbative QCD are big enough and in the right direction 
to explain the observed phenomenology at low energies. 

A very popular model for instance is that of a “freezing” coupling constant [311. Although 
certainly economic and rather successful, this approach (at least in its most naive version) runs 
into conflict with spontaneous chiral symmetry breaking as we shall show with a particular 
example in the Appendix B. 

Another alternative approach to low-energy QCD is the one which has recently appeared 
in refs. |32|]. These authors discuss an interesting attempt at describing spontaneous chiral 
symmetry breaking in QCD using a variational approach on resummed perturbation theory. 
We are curious to know how this variational approach could be related to the properties of 
large orders in the coupling constant, i.e. to renormalons, and hence to some of the results 
that will be obtained in this paper. 


The aim of this article is to study more closely and within a specific class of QCD diagrams 
the interplay between the insertion of four-fermion operators and the leading UV renormalon 
contributions; as well as the possible impact of this relationship on bridging the gap between 
QCD and the low-energy chiral effective Lagrangian. This we do by explicitly studying UV 
renormalon effects in two-point functions of colour-singlet vector currents and of colour- 
singlet pseudoscalar currents. The currents which we consider are light quark currents of 
the flavour SU(3)l x SU(3)r group. Their associated two-point functions correspond to 
physical observables in hadron physics. We are not so much interested in UV renormalons as 
a source of possible 1 /Q 2 like ambiguities in two-point function QCD sum rules, but rather 
in their possible relevance to genuinely non-perturbative effects which they may signal when 
the external Euclidean momentum Q 2 in a two-point function is taken much smaller than the 
characteristic QCD scale. With the restriction to a one-loop (3 function this means Q 2 < A|, 
where A^ is the Landau pole. Although this is not the conventional situation (wherein Q 2 is 
always assumed to be very large), we think that taking Q 2 small is unavoidable if one ever 
wants to make any contact with a low-energy effective Lagrangian for QCD, as this is the 
range of momentum which the effective Lagrangian is supposed to describe. We find that 
the ambiguities generated by the leading UV renormalon in this regime of low Q 2 hint at the 
existence of non- perturbative effective four-fermion local operators to describe low-energy 
physics, much the same as the study of the ambiguities generated by the IR renormalons hint 
at the existence of vacuum condensates of composite operators in the OPE. Moreover, these 
non- perturbative contributions turn out to be exactly like those of the extended NJL model 
mentioned earlier, with the four-fermion operators normalized by an energy scale which is 
momentum-independent, i.e. a constant. This constant which plays the role of a cutoff in 
the Euclidean momentum integrals appears to be related to the characteristic QCD scale and 
it seems natural to identify it with A x , the chiral symmetry breaking scale [33]. At the level 
of a one-loop /? function this identification implies that A x — A^. 

Most of the technical part of the paper has to do with the calculation of the contribution 
to the residue of the leading UV renormalon coming from the effective charge exchange of one 
and two-chains of gluon self-energy-like graphs in these two-point functions. The concept 
of a QCD effective charge at the one loop level is reviewed in Section where we also define 
the QCD “amputated” action which we use as a calculational framework. For the sake of 
simplicity we work in the chiral limit where the light u, d, s quark masses are neglected. In this 
limit the two-point functions of flavour non-singlet axial-vector currents and scalar currents 
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are trivially related to those of vector currents and pseudoscalar currents respectively; and 
it is therefore sufficient to study the two types of two-point functions which we do. Since, 
eventually, we are interested in a comparison with an “all-orders” analysis of QCD in the 
larg e-N c limit, we have kept track explicitly of the N c factors which appear at the various 
stages of the calculations. The calculations reported in Section [| are made in two ways. 
One method uses the Gegenbauer expansion technique of conventional Feynman diagrams; 
the other makes use of the operator product expansion technique following refs. We 

find that the dominant role played by the dimension-six four-fermion operators is indeed 
universal and comes about in the same way in the two channels which we have studied. The 
specific discussion of the calculations at low Q 2 values is done in Section [|, and it is followed 
by Section [|], which is dedicated to the conclusions and outlook. 


2 The QCD Effective Charge and Renormalon Calculus. 


In QED, the infinite subset of radiative corrections summed in the Dyson series generated by 
the one-particle-irreducible vacuum polarization self-energy function n#(A: 2 ) defines an effec¬ 
tive charge which is universal, gauge-, and scheme-independent to all orders in perturbation 
theory: 


«ef# 2 ) = £ 


4 1 _ e 2 1 


( 2 . 1 ) 


4vrl + n R (k 2 ) 47rl + ri(fc 2 )’ 

where e and II(A; 2 ) denote bare quantities. The extension of a similar effective charge con¬ 
cept to QCD is of fundamental importance for renormalon calculus, if one wants to identify 
unambiguously the infinite subset of gluon self-energy-like radiative corrections that one 
is summing in the replacement of a gluon propagator by a so-called “renormalon chain”. 
Recently, there has been substantial progress in this direction. The theoretical framework 
which has enabled this progress is the so-called pinch technique [j34-37]. The pinch technique 
is a well-defined algorithm for the rearrangement of conventional gauge-dependent one-loop 
n-point functions to construct individually gauge-independent “one-loop” n-point-like func¬ 
tions. This rearrangement of perturbation theory is based on a systematic use of the tree 
level Ward identities of the theory to cancel in Feynman amplitudes all factors of longitudinal 
four-momentum associated with gauge fields propagating in loops. In the case of QCD, the 


resulting effective charge has the form [36| 


a eff ) 


9r 


4vr 1 _ n fl (fc2) ’ 


( 2 . 2 ) 


where k ^ = —k 2 > 0 for A’ 2 -spacelike, and II/j(A: 2 ) is a gauge invariant gluon self-energy-like 
two-point function, which in the MS -scheme in particular, and at the one loop level is given 
by the expression 
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(2.3) 


The coefficient ^p-N c + of the logarithmic term in this equation is precisely the first 
coefficient (3\ of the QCD /3-function. The effective charge encodes therefore the physics of 
the /3-function; in particular the scale breaking property. Since, eventually, we are interested 
in the appearance of physical scales in QCD, it seems natural to focus our attention on the 
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properties of those diagrams of perturbation theory generated by the insertion of effective 
charge-ex changes. Of course, we do that at the simplest level of keeping only the one loop 
dependence of the /3-function. As far as one is only interested in general qualitative features 
this should not be a serious limitation. 

With the concept of an effective charge in hand it is possible to define a framework to do 
QCD renormalon calculations in the sector of light quark flavours. It is this framework which 
we next describe. 


Let T(v, a, s,p) be the full QCD generating functional of the Green’s functions of quark 
currents in the presence of external vector v, axial-vector a, scalar s, and pseudoscalar p 
matrix field sources: 


jr(v,a,s,p) 


VGu exp 


-i j dx \G^(x)G^(x)} x 


VqVqexp ^ i j d 4 x q(x)i pq(x)^j , 


(2.4) 


where q = ( u,d,s ), and p denotes the QCD Dirac operator in the presence of the external 
sources 

P = + ig s Gp,(x)\ - (x) + 75 a M (x)] + i[s(x) - i^pfx)} ; (2.5) 

Gn is the gluon gauge field colour matrix (a = 1,2,3,... 8) 


G„(I) ^ V 


( 2 . 6 ) 


and 

G$(x) = d„G^ - d v G& - gsfabcG®G® , (2.7) 

the eight gluon field strength tensor components. The factor Z in eq. (2.4) is such that 

T(0,0,0,0) = 1. 

The QCD contributions which we shall consider in the renormalon calculations presented 
below are the ones generated by an “amputated” generating functional e ir ( v,a,s ’ p ' ) defined as 
follows: 


JF(v,a,s,p) 


1 

z 


VqVq exp \i J d 4 x q(x) [i^ (d^ - i{v^ + 75o M )) + i(s - ij 5 p)] q{x) 

( r \(a) \(b) "j 

-SfiE J d * xd4 y Qixh^—qix) i&n lw (x - y ) q(y)Y—q(y) \ ■ ( 2 - 8 ) 


with iA^ w (x — y ) the Fourier transform of the renormalized gluon propagator-like function: 


k^k u 


^ { (-»» + k 2 

and d R {k 2 ) the effective charge defined in eq. (2.2) i.e. 


d R(k 2 ) -£ 


kpfkv 

k 2 


(2.9) 


9R d R(k 2 ) = -= 4vra eff (fc|). 
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( 2 . 10 ) 
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The function iA r ^ fW (x—y) results from having integrated out those gluonic interactions which, 
order by order in perturbation theory, contribute to the QCD effective charge as successive 
powers of gluon self-energy-like two-point functions evaluated at the one loop level. In terms 
of a characteristic scale, like e.g. the scale of the M 5-scheme, we have 

«eff( fc l) = —;-p— > ( 2 - n ) 

T>F lo S A\l== 

MS 


where 

and [34-gg c 2 = exp|^^:J. (2.12) 

The numerical value of the constant c in different relevant limits for N c and rif does not 
change much: c(N c = 3, n/ = 3) = 2.87 ,c(N c —► oo) = 2.76 , c(n/ —>• oo) = 2.30. We shall 
often refer to the scale 

A l = c X , (2-13) 


as the Landau pole. Notice that, with ~ 300 — 400 MeV, as determined phenomeno¬ 
logically, A l turns out to be numerically of the same size as the chiral symmetry breaking 
scale: A^ ~ A x ~ lGeV . Although c and A-^ are scheme dependent, the combination A/, 
is scheme independent. For simplicity we shall work in a subtraction scheme in which the 
corresponding c is unity, or equivalently, wherein < 7 r(/u) = gj^g(/x/c)-, and from now on we 
shall drop the subscript “R” from gR(fji). 

We consider that the framework described above is a net improvement with respect to 
previous approaches to renormalon calculus in QCD. The existence of an effective charge in 
QCD, with properties analogous to those of the QED effective charge , as recently emphasized 
in refs. [361 and 137], provides the basic feature that allows us to select a minimal class of 
well defined contributions to QCD renormalons. In the so called naive non-abelianization 
procedure (NNA) [|38|], it is suggested to perform first QED-like calculations —ignoring the 
non-abelian gluonic interactions— and replace in the results thus obtained the number of 
light quark flavours, nj, coming from QED-like vacuum polarization insertions, by 


n f 



(2.14) 


Another procedure often adopted is to consider abelian-like gluonic interactions dressed with 
the asymptotic running coupling at the scale of the virtual gluonic Euclidean momentum, 
and then advocate a “large /3\ expansion”. [] None of these prescriptions, included the “am¬ 
putated” generating functional in ( |2.8|) which we are proposing, has been justified as yet in 
terms of a well defined approximation within QCD itself. The merit however of the “ampu¬ 
tated” effective action in is that it selects a minimum set of a well defined class of QCD 
contributions; and in principle it could be improved by considering higher loop contributions 
to the gluon self-energy-like function 11^ which governs the gluon propagator-like function 
iA'pf /jIy (x — y)\ as well as by taking into consideration more and more non-local interaction 
(three-point-like, four-point-like, ...) terms. 


As an illustrative example, we dedicate the rest of this Section to the calculation of the 
contribution to the Adler function induced by the renormalon effects which result from the 

3 See e.g. ref. [^9| and references therein. 
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exchange of an effective charge chain in the “amputated” generating functional T(n, a, s,p) in 

Q. With 

i (: u^u : — : d^d :) (2-15) 

the vector-isovector quark current, and 


U^(q) = iJd 4 xe iq - x {0\T{V^(x)V u (0)}\0) 

= -(/V-?Y)n(? 2 ) (2.16) 

the associated correlation function, the Adler function is defined as the logarithmic derivative 
of n(q 2 ): ( Q 2 = — q 2 , with Q 2 > 0 for g 2 -spacelike) 

AQ 2 ) = -Q 2 f^ ■ < 2 - 17 ) 


In perturbative QCD calculations, Q 2 is supposed to be larger than the characteristic QCD 
scale, which in our case is A 2 in eq. ( 2.13| ). In terms of Feynman-like diagrams, the calculation 
in question can be represented by diagrams like the one in Fig. 1: 



Fig. 1 One of the three Feynman-like diagrams with an effective-charge chain which induce 
renormalon effects in the Adler function. 


The chain of bubbles in Fig. 1 corresponds to the replacement of the ordinary free gluon 
propagator by the full gluon propagator-like function 


k 2 + ie 


( Ws) ^ i 



kgXA 4vr a e ff( fe l) 
k 2 ) — ie 


-t- kphy (- ig s ) 2 

k 2 k 2 + ie 


(2.18) 


where kE denotes the Euclidean virtual momentum carried by the chain; and £ is the same 
covariant gauge parameter as in eq. (2.S). This replacement is in fact the net effect, in 
momentum space, of the interaction term in eq. (2.8) when evaluated to its lowest non trivial 
order. 

An interesting function to consider then is the functional derivative 


SIl{Q 2 ) 

Hd#!) ’ 


(2.19) 
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Up to overall normalization factors, it corresponds to the forward elastic scattering amplitude 
of an off-shell vector-isovector quark current off an off-shell gluon evaluated at the one loop 
level as illustrated in Fig. 2. 


k E k E 



Fig. 2 One of the three Feynman diagrams corresponding to the forward scattering ampli¬ 
tude of an off-shell vector-isovector current, carrying Euclidean momentum Q, off an off-shell 
gluon carrying Euclidean momentum kE- 


Once this function is known, the full calculation of renormalon effects induced by one effective 
charge chain exchange is simply given by the integral 


n(Q 2 ) 



SU{Q 2 ) 
6a eS ( fc i) 


®eff ^ e ) ■ 


( 2 . 20 ) 


It turns out that the kernel given by the functional derivative above can be extracted from 
early papers on QED by Baker and Johnson |40|]. It can be easily adapted to QCD by 
including the appropriate colour factors and it has the following form ( Cf 


2N C 


jn(Q 2 ) 


N c 

32t r 3 


Cf k^ x 


fF S (w) for - Q 2 ’ 
for k E ^ Q 2 ’ 


( 2 . 21 ) 


clearly illustrating the UV IR symmetry of the kernel. Notice that renormalization does 
not affect the calculation of this kernel. At this level QCD with massless quarks still has its 
full conformal invariance properties, which is reflected in the symmetries of the kernel above. 
The function H(z) has the simple integral representation 


■(z) 




1 

i + y 


log y. 


( 2 . 22 ) 


In the interval 0 < z < 1 it is a monotonically decreasing function from S(0) = 1 to 5(1) = 
^ — rf7T 2 = 0.2802..., and it has the Taylor expansion: 



l2 log " 


13 

144 


4 /1 

F(z) = l + -z ( -i°g« 


+ ••• ■ 


(2.23) 








For the purpose of renormalon calculations, it is natural to separate the integration region 
in eq. ( 2.20 ) into an infrared-dominated region where 0 < k E < Q 2 , and an ultraviolet- 
dominated region where Q 2 < k E < oo. The integral in the infrared region is in fact ill- 
defined because the function a e f[(k E ) blows up at the Landau pole where k E = A 2 . [] As a 
consequence of this, there is an obstruction as a matter of principle to reconstructing the full 
non- perturbative answer from resummed perturbation theory. This intuitive argument can 
be refined and made more systematic with the use of Borel resummation techniques Q . Let 
us discuss the most salient features of the integral in eq. (2.20): 


The limit a c ff{k E ) = a(/t 2 ). 

In this limit we find the well known perturbative result corresponding to the one gluon 
exchange contribution to the Adler function, namely 




(2.24) 


Notice that in this limit the integral in the region from 0 < k E < Q 2 contributes a 
constant to n(Q 2 ) and hence nothing to the Adler function; i.e. it is renormalized 
away. One-gluon exchange with no gluon self-energy-like correction (no bubble) in it 


t? Q 


0 = 1 in the 


is not sensitive to the infrared region. Only the first term ^ . , 2 

\ k e 

expansion of the S function in the UV region, where k E > Q 2 , contributes to the Adler 
function in this limit. 

Infrared renormalons. 

Once we start summing bubbles corresponding to the expansion of a e g(k E ) in powers of 
gluon self-energy-like insertions, we find that the leading infrared renormalon is induced 

by the first term, E I —>■ 0 1 = 1, in the expansion of the E-function in the IR region. 

Making the change of variables 


w/2 = -b 0 a(Q 2 ) log k\/Q 2 
in eq. ( |2.20D , one obtains 


where 6 q = — /3i/2it 


n (Q 2 


lIR 


^ Cf 1 


dw e b o a (d 2 ) 



167r 2 27 t6o Jo 

which leads to the following contribution to the Adler function 


(2.25) 


(2.26) 


N, 


1 


-* W)l ® = i6F CF W„ dwe 



(2.27) 


This expression is already in the form of a Borel transform. An expansion around w = 0 
would generate the characteristic n! behaviour of the perturbative expansion in a(/r 2 ) 
which we mentioned in the Introduction. As already discussed by other authors &0 

4 This is in fact only true for a one-loop f3 function. In the case of a two-loop f3 function, the relevant pole 
is not at tif [ll], but the physics is pretty much the same. 
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|42|j , we find that the leading IR renormalon contribution to the Adler function appears 
as a pole in the Borel plane at w = 2. There is no term in the IR expansion of the 
k 2 \ 

I function which leads to a pole at w = 1. From the previous change of variables 

one also sees that low values of w are associated with momenta around the large scale 
Q 2 , where perturbation theory is expected to give a good description of the dynamics. 
However, as w goes up (and it has to go all the way up to infinity) one enters deeper 
and deeper into the IR region, where perturbation theory must fail or else QCD would 
not describe the spectrum of hadronic bound states that are observed. The singularity 
at w = 2 exhibits this in its crudest form. More precisely, this singularity implies an 
ambiguity in the perturbative evaluation of the Adler function. Therefore, the analytic 
continuation in w in eq. (2.27), or equivalently the resummation of perturbation theory 
into the effective charge Oi^{k 2 E ) of eq. ( 2.2C| ), that one has tried in order to obtain the 
full solution has failed. The ambiguity is encoded in the unavoidable prescription to 
skip the pole; and is of the form ~ e -i/&o<*(Q) . 


(L4(Q 2 )Iir 


N c C F 


167T 2 —fli 


Al \ 4 

Q) 


(2.28) 


i.e., an ambiguity which has the same -X% pattern as the gluon condensate contribution 
that appears in the OPE evaluation of the Adler function [|T|, 131. Since the Adler 


function must be an unambiguous physical observable, there must exist another con- 

In other words, (Borel) resummed perturbation theory 


tribution that cancels 
requires the presence of the gluon condensate (which is also ambiguous for the same 
reason) to combine with the result of eq. ( 2.27 ) and yield a final well-defined answer ||. 
This is an example of how all-orders perturbation theory is capable of “hinting” at 
non-perturbative dynamics. 


Following analogous steps, the insertion of higher powers in k^/Q 2 from the IR expan¬ 
sion of the E-function in the 0 < k% < Q 2 integrand produces contributions to higher 
singularities in the Borel plane corresponding to higher order IR renormalons located at 
w = 3, w = 4, and so on. However, the higher power terms of the E-function only gen¬ 
erate partial contributions to the higher order IR renormalons. For example, for the IR 
renormalon located at w = 3, there will also be further contributions originating in the 
three-point-like non-local term which will appear as a correction to the “amputated” 
generating functional in ( t2.8|) . This new term will generate diagrams like e.g. the one 
shown in Fig. 3 which we are not taking into account here. 
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Fig. 3 Example of a Feynman diagram generated by the three gluon-like vertex. 


The total contribution to higher order IR renormalons is associated with higher dimen¬ 
sional condensates with a larger power in 1/Q 2 than the one in ( |2.28| ), and in this sense 
they are non-leading. [] 

Ultraviolet renormalons. 


The integral in eq. ( |2.20| ), in contrast to the behaviour in the infrared region discussed 
above, is well behaved in the ultraviolet region where Q 2 < k E < oo, provided that 
Q 2 » A 2 . Nonetheless it also has associated with it a Borel representation, the only 
difference being that now the poles will occur at negative values of the Borel variable 
w. As in the case of infrared renormalons, the further away the pole is from the origin 
w = 0 the smaller will be the n-th coefficient accompanying a(Q 2 ) n in the perturbative 
series. This gives rise to a hierarchy of ultraviolet renormalons. The leading one, within 
the exchange of one power of the effective charge , is the one induced by the O(z )-term 


in the UV expansion (k E > Q 2 ) of the S-function in eq. ( 2.23 ). With the change of 
variables w = boa(Q 2 ) log k 2 E /Q 2 one obtains 


n(Q 2 )| uv = - 


N c 

167T 2 


C F i 


1 


(1 + w)‘ 


9 27T&0 Jo 

11 1 

+ 


dwe b o a U 2 ) 


6 1 + w 


(2.29) 


where the leading term is the one induced by the term | ^^-ilogp-^ 

eq. (2.23). The corresponding contribution to the Adler function is then 


m 


N c „ 4 


- 4(0)| rv = T^ Cf 9 


1 


2 nb 0 Jo 

1 5 

+ 


dw e b o a (Q ) 


t} 


(l + w>) 2 ' 6(1 +re) 6 J ^ ^ 

where we have also scaled p 2 at Q 2 . The result of the leading term ~ agrees with 

the one obtained in refs. |^O,21,|0] using other methods. One sees that the singularity 
is at w = —1, i.e. outside the integration range. Therefore, unlike the case of IR 
renormalons, one cannot argue now that the singularity is an obstruction to do the 
integral. This is also clear from the above change of variables since for 0 < w < oo one 
never leaves the region of very large momentum where perturbation theory is supposed 
to be a faithful description of the dynamics. 

Contributions to higher-order UV renormalons can be obtained from the insertion of 
successive higher power terms in the UV expansion of the E-function in the integrand 


in eq. ( 2.20 ) for Q 2 < k 2 E < oo. The next Section is fully dedicated to the study of UV 
renormalons in two-point functions. 


5 In fact, in the particular case of the IR renormalon located at w = 3, and with no extra hard a(Q 2 ) 
correction, we expect a cancellation of the various partial contributions to the residue of the pole at w = 3, in 
accordance with the fact that there is no three gluon condensate contribution in the OPE to the Adler function 
as proved in refs. 0Q. 
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3 UV Renormalon Contributions to Two—Point Functions. 


We shall discuss in this Section the calculations of the contributions to two-point functions 
of quark current operators induced by the leading UV renormalon generated by the exchange 
of one and two QCD effective charge chains. The two-point functions which we shall consider 
are the Adler function already discussed in the previous Section, and the two-point function 
associated with the divergence of the axial-current: 

d^Affx) = (m d + m u ) : d{x)i^u{x) :, (3.1) 

where the overall quark mass factor will only be kept so as to make the pseudoscalar current 
density renormalization invariant. All the calculations presented here have been made in the 
chiral limit where the masses of the light u, d, and s quarks are set to zero. The two-point 
function associated with the divergence of the axial current is 

T 5 (g 2 ) = i J d A xe iq ' x { 0 | T {d»A^x)d u At (0)} | 0). (3.2) 

The equivalent of the Adler function here is the second derivative of Ts(g 2 ), 

Q 2 8 2 T 5 (q 2 ) 

~~ (m u + m d ) 2 ( dQ 2 ) 2 

This second derivative is, like the Adler function in eq. @, independent of external sub¬ 
tractions. 


V{Q 2 


3.1 The Vector Two—Point Function. 


As mentioned in Section [| we have carried out the UV renormalon calculations using two 
different methods. One is the operator product expansion (OPE) technique, which we discuss 
below, following rather closely the QED work of Vainshtein and Zakharov [20, [21]. The other 
method uses ordinary Feynman diagrams directly and the Gegenbauer polynomial expansion 
of Euclidean propagators. For the vector two-point function, we shall discuss this second 
method only in connection with the exchange of two effective charge chains. 


3.1.1 OPE Calculation: One Chain. 


Let us start with the contribution from the UV renormalon generated by the exchange of 
one effective charge chain, to which we shall refer for short as the contribution of a “one- 
renormalon chain” or even “one-chain”. The contribution of the “one-renormalon chain” to 
the external vector field vacuum polarization tensor of eq. ( [2.16| ) which we have discussed in 
Section |2] can also be written as follows: 


n ^(q) v»v v 


-i f d A k 47ra ff (fc|) T 

T (Sp < v 1 F 1 ' 


(3.4) 


where the external vector field v is the same as the one appearing in the Dirac operator in 
eq. (2.5) and T is the time-ordered operator 


T = ik 2 / d A xe 


ik-x 


T 


: q {x) ■ ■ <?(0)7My9( 0 ) : 


(3.5) 
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Following |2^,^T| we now perform the operator product expansion on the operator T by using 
the Schwinger background field technique, explained e.g. in ref. [[SJ. The first nonvanishing 
contributions comes from the dimension-six operators: 


t = W2 


fleff^j) 


2 Nr. 


D a G afl (0) + C F D a F a J0) 


9(0), 


(3.6) 


where F /lll (x) is the field strength tensor associated with the external vector source. With this 
expression for the operator T one can easily compute the integral in (30). [Recall eq. ( 2.16 ), 
and notice that at this order only the F afl term in (3.6) contributes.] This results in the 
expression 


n(g 2 ) = - 


N c 2 
16-7T 2 9 


Q 


IQ 2 


dk 2 E 

k% 


a eff(^ 


7 r 


log(A:|/Q 2 ) , 


(3.7) 


where the momentum in the argument of the logarithm has been cut off at the scale k E that 
appears in the OPE in eq. (|3.6|). The steps of the calculation are illustrated in Fig. 4 below: 




Fig. 4 The leading UV renormalon contribution from a one effective charge exchange, 
evaluated with the operator product expansion technique, is the one generated by the dimension 
six vertex like operators simulated by the form factor in the lower figure. 


The UV renormalon contribution to the Adler function in can now be readily 

obtained. Using the change of variables w = boa(Q 2 ) log(k E /Q 2 ) results in 


AQ 2 ) 1^ 


"" c F 4 


1 


dw 


I 67 r 2 9 27r&o Jo (1 + w ) 2 

which, if expanded in perturbation theory, yields 


g b 0 a(Q 2 ) 


A(Q 


2\ 1 one chain 

luv 


- E 

Large n 


Nr 


;C F - 


4 1 


167T 2 9 27t6q 


(- 6 0 a(Q 2 )) n\ 


(3.8) 


(3.9) 


Notice that the result in (3.8) is in agreement with the one we found in eq. ( 2.3C| ) for the 
leading behaviour of the Borel transform around the singularity w = — 1. Our result coincides 
also with the one obtained in refs. [|(], 211 in the abelian limit where ^-> 1 . 
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3.1.2 OPE Calculation: Two Chains. 


In an analogous way one can deal with the more complicated contributions coming from “two 
renormalon chains”, as illustrated in Fig. 5. Two different new structures emerge from the 
box-like diagrams: 


T 


Box 


t-Box , /-t-Box 
1 1 + 1 2 


r_/„x n A a^ a 

#hV7 a T T^ 


(3.10) 





Fig. 5 The leading UV renormalon contribution from the exchange of two effective charges. 
The evaluation with the operator product expansion technique, generates a dimension six four- 
fermion operator simulated by the box form factor in the lower figure. 


Using the operator identity 

(7 m 7 A 7“) 0(7a7A7/x) = 10 Y l 0 7m - 6 7^75 0 7m75 , 
one can immediately see that the abelian limit yields the result 

6 g 2 a{k 2 F ) 

effV EJ 9(0)7^59(0)9(0)7^759(0); 


X; 


Box 


abelian 


k 2 


(3.11) 


(3.12) 


i.e. there is a cancellation of the 7 ^ structure. Equation ( |3.12| ) was also obtained by the 
authors of refs. [20, 21 ]. However, in the large--/V c limit, the previous cancellation does not 
take place. Instead one sees that 72 ® ox is proportional to N c whereas 7)® ox goes like 1 /N c , 
i.e. it is suppressed by two powers of N c relative to XP° X . Consequently, at large-lV c , 


T 


Box 


~ To 


Box 


-jVc 9 eS (k 


4 k 2 


10 ?(0)7 /i y«(0)g(0)7 M yg(0) 


-6 9(0)7^75^-9(0)9(0)7^75 ^-9(0) 


(3.13) 
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Upon comparing with the one-renormalon chain contribution of eq. (3.6) we see that the 
use of the equations of motion for the gluon field in that expression also produces a four- 
fermion operator of the form ( ]3. 13 ) 2(], 21] [] which is however suppressed relative to the one 
in eq. ( |3.13 ) by a factor 1/./V 2 and hence we neglect it. Therefore, we find that large—Ay selects 
as the only operator contributing at the “two-renormalon-chain” level the one in eq. 

If we now Fierz the operator TP 0X , and use the identity 


a ~ b = 


a — b 


E I q b L I 2 +(L - R) 

a,b 


- 2(a + b)E(^)fe£), (3-14) 

a,b 


we find that the result is almost like the four-fermion operators which appear in the ENJL 
model of large- N c QCD [28,^]. There is however a very important difference, namely that 
the dimension-six T® ox operator which appears here is suppressed by the momentum scale 
k E and not by a constant scale like A 2 . We shall come back to this important issue in the 
next Section. Here we shall limit ourselves to a numerical observation concerning the relative 
size of the scalar-like coupling Gs ~ 2(a + b) and vector-like coupling Gy ~ (a — b)/2 which 
appear in eq. (|3.14|) . Plugging in the values a = 10, and b = 6 obtained in eq. ( 3.13|) results 
in a ratio Gs/Gy = 16. For comparison we recall that a cut-off one-gluon exchange yields 
b = 0 ,Gs/Gy = 4, and the abelian version of refs. 20|, 21], which has a = 4/3, b = —3, 
would yield Gs/Gy = —20/13, a negative value !. The phenomenologically favoured value of 
ref. |28] is Gs/Gy ~ 1. 

Let us come back now to the contribution of “two renormalon chains” to the Adler func¬ 
tion. The relevant diagram is the one depicted in Fig. 5, where the black square stands for 


the insertion of the operator ( 3.13 ). As already noticed in [20, 211 a subtlety arises at this 
point, for if one wishes to keep the relations (3H) and ([ho]), in order to avoid double counting 
one has to multiply the right hand side of eq. (EO) by a factor 1/2. After Fierz-ing the result 


in ( [3.13 ), only the vector part contributes of course, and one obtains (recall eq. ([2.16 ) and 
the fact that we are working in the large- N c limit) 


nfa 2 ) = Nc 




16-7T 2 V 2 27 tJ V 9 J Jq 2 kg 


/‘OO si 

/ ahj E 


a eS( k E )f l °S 2 {k 2 E /Q 2 


(3.15) 


for the external vector vacuum polarization function, and 


A(Q 


2 \ |two chains. 
IUV 


N r 


N c 1 y 2 


16vr 2 V 2 27T&Q / 9 


0 

2 ' 2 \ 


dw e b o a (Q 2 ) 


for the Adler function, where w = boa(Q 2 ) log(A ; E /Q 2 ). 
At large orders one finds 


w 


(1 + w) 3 


(3.16) 


A(Q 2 ) chains ^ 


N c 1 \ 2 1 


N r , 


167T 2 V 2 27T&0 


) g E (- fo oa(Q 5 


Large n 


(n + 1)! 


(3.17) 


It turns out that, both expressions ( |3.9|) and ( [3.17 ) are comparable in the large-W limit and 
of the same order as the one-loop “parton” graph. However, ([3.17) is leading at large orders 


6 i.e. a penguin- like diagram. 


15 
























































of perturbation theory, i.e. for large n. In this sense one can view the “two renormalon 
chain” contribution as a selection of a subset of the whole set of diagrams that are leading 
at large-lV c . This simplification may be welcome in the sense that, so far, even the leading 
contribution at large-lV c has proven itself to be already too difficult to deal with in QCD. We 
shall comment on this again in the next Section. 


3.1.3 Two Chains with the Gegenbauer Expansion Technique. 

We want to present another technique for calculating the “two renormalon chain” contribution 
to the Adler function which uses directly the Feynman diagrams generated by the “ampu¬ 
tated” effective action in (2.8). Further limitation to leading contributions in the large-lV c 
limit, and the fact that we are only interested in the leading UV renormalon contribution 
restricts the class of possible diagrams with two chain exchanges to the one in Fig. 6, with a 
Feynman gauge-like coupling i.e., in this diagram only the g^ v term in (2.9) is operational. 



q 


Fig. 6 Diagram representing the exchange of two effective charges. A convenient routing 
of momenta to evaluate the contribution to the leading UV renormalon using the Gegenbauer 
expansion technique is indicated in the figure. 


With the routing of momenta indicated in Fig. 6, the contribution to the vector two-point 
function reads as follows: 


m^(q) = N c 


Nc 1 

2 2 


d 4 k\ f d 4 k 2 


d 4 k 3 


(4vr) 2 


dz\e 


-Z 1 


(2vr) 
1 


(2tt) 

X3 

dz 2 e 


■x 


(2vr) 4 ' 

g aa '(h 2 ) b0Z1 9w>(h 2 ) b0Z2 

(k 2 — k\) 2 A+ b o z i) (^ 3 _ ^ 2 ) 2 ( 1+fe ° 22 ) 


tr 


Y fti'y 01 fo'yP hY{h~ 4)l a \h~ 


klk%k%(ki - q) 2 (k 2 - q) 2 (k 3 - q) 2 
where we have used the identity 


(3.18) 



dze 


i 



(3.19) 
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In order to calculate this integral, it is convenient to first introduce Dirichlet variables Q: 

W 2 


m = b 0 zi 


w 2 = b 0 z 2 


w = W\ + w 2 , 


V = 


w 


(3.20) 


and expand the denominators which depend on differences of momenta in Gegenbauer poly¬ 
nomials rp 

u 1 -1 

(3.21) 


{k p) 2X (fc| - 2k E PECOs9 +p E ) 2X 


(3.22) 


n=0 


where the C^(z), [C^(z) = C n (z)] are Gegenbauer polynomials with the orthogonality prop¬ 
erty: 


dn k C*(a • k)C^ n (b • k) = S n , r 


A 


A + n 


C^d • b) 


(3.23) 


with Cn( 1) = "Tnrpxp and dkl k the solid angular element resulting from 

r 2tt ( 2_e / 2 ) 

J d (4 ~ e) k - i fe) 3 “ £ . 


(3.24) 


The trace in the numerator of eq. ( [3.18| ) also has to be expressed in terms of powers of the 
Euclidean momenta k\E, k 2 E, k^E, Q = QE\ and powers of cos 0*0 = cos (ki ■ q), i = 1, 2,3; and 
cos 6ij = cos {ki ■ kj ). The Dirac trace results in at most three powers of cos 0’s; and therefore, 
from the orthogonality properties of the Gegenbauer polynomial, only a few terms give non¬ 
zero contributions to the solid angle integrals / dQ kl dfl k2 dfl k3 . With further restriction to 
those terms which can contribute to the leading UV renormalon, one finds that only terms 
with two powers of ^ and terms with no powers of in the Dirac trace are relevant, with 
contributions proportional to Q^k^E and — | Q^k^E respectively. This results in a contribution 
to II(g 2 ) from the leading UV renormalon generated by the exchange of two chains: 


n (q 2 ) = 



1 

2 27T&0/ Jo 
„!.2 ,, 2 


dw w e 6 o“OD x 


rk 2 
I h '2E 


dk 


IE 


IQ 2 


k 2 

He 


p u2 
I hj 2E 


’Q 2 


^3 E 


91 

h 2 


J-X 

2 nbo) 


, _ i / n 2 \~ w r i i 

dwwe Wb o a (X 2 ) — / dxx w log 2 — . (3.25) 

\ H J Jo x 


The corresponding contribution to the Adler function, scaled to p 2 = Q 2 , is then 


^(Q 2 )uv chains 


N c (N c 1 \ 
167T 2 V 2 27r6o / 




dw e Wb o a (Q 2 ) 


(1 + wY 


(3.26) 


7 The change to the Dirichlet variables w and v automatically implements the convolution property of 
Laplace transforms. 

8 Useful properties of Gegenbauer polynomials in connection with Feynman diagram calculations can be 
found, e.g., in ref. jd^] . 
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in agreement with the OPE result in eq. (3.16). 


The final result of the leading UV renormalon contribution from the exchange of one and 
two powers of the QCD effective charge to the Adler function A(Q 2 ), i.e., the sum of one and 
two renormalon chains, is then given by the expression: 


A(Q 2 ) 


Nr Nr 14 


16vr 2 2 2v rb 0 9 J 0 


dwe U 6 o“(<3 2 ) 


1 


Nr 1 1 


1 


(1 + u;) 2 2 27r6o 2 (1 + w) 3 


(3.27) 


where the first term in the second line corresponds to the contribution from the one renormalon 
chain, and the second term to the contribution from the two renormalon chain. 


3.2 The Pseudoscalar Two-Point Function 

We shall next discuss the contribution of the leading UV renormalon to the two-point function 
^(g 2 ) defined by eqs. (|T2|) and 0 , calculated also at the level of one and two chains. 
Here we shall first describe the calculations we have done using the Gegenbauer expansion 
technique. 

3.2.1 One Chain with the Gegenbauer Expansion Technique. 

There are three possible diagrams with a one renormalon chain contributing to the two-point 
function Ts(g 2 ). Much the same as in the case of the vector two-point function, if we restrict 
ourselves to the contribution from the leading UV renormalon, it is then sufficient to consider 
the contribution from the one renormalon chain exchange in Fig. 7 with a Feynman gauge-like 
coupling of the renormalon chain. 



Fig. 7 Diagram representing the exchange of one effective charge. A convenient routing 
of momenta to evaluate the contribution to the leading UV renormalon using the Gegenbauer 
expansion technique is the one indicated in the figure. 


With the routing of momenta indicated in Fig 7, the contribution to ^(g 2 ) in large-_/V c is 
then as follows: 


^5 (q 2 ) 


—i(m u + m d ) 2 N c 


N r 


d 4 k i 

W) 


dfik -2 -i (~ig e ff [~(fci ~ fci) 2 ]) 2 
( 27 t 4 ) (hi - k 2 ) 2 + ie 
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tr[i75 hi 11 failbifa- 


[k\ + ie] [k% + ie] (hi - q ) 2 + ie (k± - q) 2 + ?'e 

Using the identity 

~i (~ig e ff [~(fci - fc 2 ) 2 ]) 2 _ .4vr 


(3.28) 


(A:i - k 2 ) 2 + ie 


= — i— tt / (fee 
/U 2 Jo 




M 


l+ 6 0 2 


(&1E — fee) 2 


(3.29) 


and expanding the denominators which depend on differences of momenta in Gegenbauer 
polynomials as in subsection 3.2.1, one obtains the expression 


T , 2\ (m u + m d ) 2 N c N c 1 
[Q ) = - 




167r 2 2 27 t6o Jo 


dwe Wb o a (k 2 ) 


1 +tU 


/Q : 


°° dk\ E ( Q 2 
2 k\ E \ k 2 


1+W 


IE t 


Y c ni (cos e 


10 , 




Q 2 

ni 


;* §p / <«*, / ((g)' 


\k 2E 


i 

n 2 


^C n2 (cOS 6*20 ) 

n 2 

16 [fc 2 E fcf E - k\ E k 2 EQ cos 0 2 o - k% E kiEQ cos 0i O + k±Ek 2 EQ 2 cos 0 1O cos 0 2 o , (3.30) 

where to = b$z and cos 0*0 = cos (hi ■ q) for i = 1,2 and cos$i 2 = cos(fci ■ k 2 ). 

The leading UV renormalon contribution in eq. ( 3.30 ) comes solely from the l = 0 term 
and the k\Ek 2 EQ 2 cos 0\o cos $20 term resulting from the Dirac trace. The angular integrals 
can then be trivially done, with the result 




(m u + m d ) 2 N c N c 1 


/T 


167T 2 2 2irbo 


(^) 


l-\-W 


°° dk\ E ( Q 2 


1 +W 


k 2 


We) 4QM ° g # ; 


Iq 2 k\ 


IE 


(3.31) 


and therefore 


*5(« 2 )Iiw ch “” = (™« + «y) 2 M 2 K Ne 1 


16-7T 2 2 27t6o Jo 


> —y-y / Q 2 

dw e ) — 

V 


1 — W 


(i + wy 


(3.32) 

Taking two derivatives with respect to Q 2 and scaling p; 2 at /r 2 = Q 2 we get the corresponding 
contribution to the pseudoscalar-pseudoscalar V(q 2 ) correlation function: 


tVCv chain 


yye± j_ r 

167T 2 2 27T&0 Jo 


dw e b o a (Q~) 


(1 + w) 


2 ' 


(3.33) 


3.2.2 Two Chains with the Gegenbauer Expansion Technique. 


As in the case of the two-chain evaluation of the Adler function in subsection 3.1.3| , the rele¬ 
vant diagram is the one in Fig. 6 with now external pseudoscalar sources, and with a Feynman 
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gauge-like coupling in the renormalon chains. With the routing of momenta indicated in Fig. 6 
the contribution to T 5 (q 2 ) can be readily written as follows 


^ 5 Or) = i( m u + "o) 2 N c 


NA 2 . n 2 


7 r 7 r 

1 — 7 , i 


r 00 Jh .2 
' a ^2E 


k 2 E dkf E f k 2 E dk 2 


3E 


(27r) 4 ( 2 vr ) 4 ( 2 vr ) 4 Jq 2 k 2 E Jq 2 k 2 E Jq 2 k 2 E 


x / / dflfcg tr [*75 ^ 17 ^ ^ 27 " 4 )lYh~ 


/4vr\ 2 f°° J -«-irry r° , 

J J dw±e fc 0 “(M J ) J dw 2 e b o a 


X ^ ni Cng (COS 6>io) 
.2\lDl+U>2 


Q \ ni v ^ r a \ f 

— I S„ 2 (7 ra2 (cos 6 * 20 ) I- 1 


ri2 


&ie/ 


V k 2E J 


Crt3 (COS 030 ) 


7_Q_ 

\fc 3 E 


n 3 


(|f) ‘ E * c £ H "(“-> (If) 3 ■ < 3 - 34 > 


The configurations with & 2 E < Aq# in the Gegenbauer expansion have been dropped because 
they do not contribute to the leading UV renormalon. The evaluation of the trace gives the 
rather simple result 


tr [...] = -64 [k\ E - ki E Q cos 0i O ) (k 2E ~ k 2E Q cos 0 2 o) (j4e ~ h eQ cos 0 3O ) • (3.35) 

To proceed further we introduce Dirichlet variables, as in eqs. ( 3.20| ). The integral over 
the variable v can be trivially done; and the contribution to the leading UV renormalon comes 
only from the terms with l\ = l 2 = n 2 = 0, and n\ = 77.3 = 1. The integral then simplifies a 
lot, with the result 


1 


(m u + m d ) 2 


^5 (q 2 ) = 


N c 

I 67 r 2 


Nc 

2 


4 Q 2 


r°° dk 2 E ( Q 


(2 irb 0 ) 2 

2 \ 1+™ 


r°° 1 

dw w e b o a (k 2 ) x 


>Q 

N, 


2 Ye 


. Ye , 


log 


k 2 ( 11 2 
2 k 2E I h 


167T 2 


N r 


4^2 


( 2 vr 6 0 ) 5 


Q 2 VQ 2 / 

1 1 

dwe b o«7 2 ) 


2 u> 


(1 + w) 3 


f) < 3 - 36 > 


Taking two derivatives with respect to Q 2 and scaling /i 2 at /r 2 = Q 2 we get as a final result 
for the two-chain contribution to the pseudoscalar-pseudoscalar V(q 2 ) correlation function 
the expression: 


7V)lw chains 


N c (N c 1 

16-7T 2 \ 2 27t6q 


(-16) 


f°° _ w _L_ 

dwe b o a (Q ) 


(1 + ie ) 3 


(3.37) 


3.2.3 The One and Two Chain Calculations with the OPE Technique. 


The calculation of the pseudoscalar two-point function with the OPE technique is done along 
the same lines as the one described earlier for the vector two-point function. Recalling the 
definitions (2.5), (3.1), (|3.2|) , (|3.3| ) one finds that 


*5 (q 2 ) 


1 f d 4 k 47TO ff (fc|) T 

2 v (p| F lri 


(3.38) 
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where T is the same as in eq. (|h^) , and | p) is the pseudoscalar “particle” annihilated by the 
external field p(x) of eq. (2.5). Contrary to what happened for the vector two-point function, 
now at the level of the “one renormalon chain” one does find a priori a dimension-four 
operator in the expansion of the operator T which, in the large-lV c limit, reads 


TdimA = 3 N c 9 ( 0 ) 2759 ( 0 ) p( 0 ); 
as well as some new dimension-six operators 

%im. 6 = pr 9(0)*759(0) 2 d 2 p( 0) + 4 p(0 ) 3 


(3.39) 


(3.40) 


besides, of course, the term D a G afl of eq. (3^) which, being subleading in N c , we shall not 
consider any further. Dimension four operators do not contribute to V(Q 2 ) | 2 (], 21 ], As for 
the contribution of the dimension six operators of eq. ( 3.40|) , the relevant diagram is the one 
analogous to the vertex form factor depicted in Fig. 4 and yields 


V(Q 2 ) chain 


4N? 


r°° dk 2 E 


Wf Q U it a * s(ki) los(ki/Q 


which, with the change of variable w = boa(Q 2 )log(k E /Q 2 ), reads 


(3.41) 


7V)luv Chain 


167 T 2 2 27 T &0 Jo 


dw e 6 o“(m 2 ) 


(1 + w) 2 ’ 


(3.42) 


in agreement with the result ( 3.33 ) obtained with the Gegenbauer expansion technique. 

When considering two renormalon chains, the four-fermion operators of eq. ( [3.13| ) also 
contribute to V(Q 2 ) through diagrams analogous to the ones in Fig. 5, with the result 


or with the above change of variables 


r™ dk% 


Q 2 


1 2 


“efffe) lo g {^e/Q ) 


(3.43) 


V{Q 2 ) chains 


N r 


N c 1 


167T 2 V 2 2nbo 


w 


16 / dw e b o a (Q ) 

Jo (1 + w ) 3 


(3.44) 


also in agreement with the result obtained in eq. (3.37) 


The final result of the leading UV renormalon contribution from the exchange of one and 
two powers of the QCD effective charge to the function V(Q 2 ); i.e., the sum of one and two 
renormalon chains, is then given by the expression: 


V(Q 2 ) 


N c N c 1 

167T 2 2 27T&0 


f 00 w 1 
dwe b o a (Q 2 ) 


{ (1 + w) 2 


N < 1 2 1 1 

2 2irbo (1 + w) 3 J 


(3.45) 


where the first term in the second line corresponds to the contribution from the one renormalon 
chain, and the second term to the contribution from the two renormalon chains. 
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3.2.4 Comments on the two—point function calculations. 

We shall conclude this Section with various comments based on the previous calculations. 


• Large-N c equivalences. 

Recall that 60 = Tjy, and /3± = -^N c + \ n f- Since we are working to leading order in 
the l/JV c - expansion, we have the following equivalences: 


1 c 1 N c 3 

27t6q F 27t6o 2 11 


(3.46) 


We keep the second algebraic form explicitly in most of the results because it shows 
better the origin of the various factors and the fact that we are neglecting non-leading 
contributions. 


• Large order behaviour. 

The large order behaviour in perturbation theory of the Adler function and the pseu¬ 


doscalar correlation function V(Q 2 ) can be read off from the results in eqs. ( |3.27 ) and 
( 3.45 ) above. They are: 


and 


-4(<3 2 )l uv ~ ^2 (ha.(Q 2 )) n . (3.47) 

p«3 2) luv ~ 1^2 IT* -11 ”" 1 "' (‘ -('»“>( < 3 2 ))" • < 3 - 48 ) 

As already discussed above, these results show that, asymptotically, the effect of the two 
UV renormalon chains is leading with respect to the one UV renormalon chain. Notice 
also that the overall coefficient of the pseudoscalar T’(Q 2 )|pjy-function is an order of 
magnitude larger than the corresponding one in the Adler A(Q 2 )|pj-y -function. Taking 
the results (3.47, jOS ) at face value, one would infer that the dominance of the two 
renormalon chain in the pseudoscalar correlation function happens at earlier n values 
than in the vector correlation function. 


• The expansion in terms of local operators. 

As we have seen from the previous calculations, one and two chains of bubbles generate 
at a large Euclidean scale a set of local operators (on the scale 1 /&#) of increasing 
dimensionality, i.e., 

r = <w>. (3.49) 

where n* = dimension of the operator 0j(O). In this equation Ci(kE ) are the corre¬ 
sponding Wilson coefficients to be obtained as a power series in oi Q ^{k 2 E ). The crucial 
observation, already made in refs. | 2 (],[ 2 l|], is that the physics starts with operators Oi 
of dimension not smaller than six. As explained in | 2 l|] operators of dimension smaller 
than six have their contribution to the two-point functions buried in the infinite renor¬ 
malization constant already present at the parton level and that drops out once enough 
derivatives with respect to the external momenta are taken in the corresponding two- 
point function. 
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Furthermore the expansion of eq. (3.49) is meaningful, in the sense that contributions 
from operators Oi of higher dimension give rise to singularities in the Borel plane fur¬ 
ther away from the origin, i.e. their n —th coefficient in the asymptotic perturbative 
expansion is more and more suppressed. 

• Diagrams with increasing complexity. 

As one considers contributions coming from diagrams with an increasing complexity, 


i.e. more and more chains, the different Wilson coefficients in ( 3.49[ ) have more 

and more powers of a c g(k E ). One can easily see that integrals like 


Ip — 


Q 2 J Q2 [« e f#i)] p lo £ 2 {k 2 /Q 2 


(3.50) 


which will appear in the two-point functions ( |3 .15 ) and (3.43), give a contribution to 
the leading term in the perturbative asymptotic expansion that is independent of p and, 
consequently, powers of are not suppressed. To see this just use 


[«ef= 


HQ 2 )} 


2 11P 


(1 + 6 0 a(Q 2 ) \og{k 2 /Q 2 )) p 

-fr 0 a(Q 2 )log(fc 2 /Q 2 ))” r(re + p) 

r(p) 


oo 

MQ 2 )F E 

n=0 


n\ 


and insert it in eq. (3.5C). After the k integration one finds 


i °° /_ l, \n 

1 = w-T E { -ff-a(Q 2 )" +r T(n + pMn + 3) . 

v* / n .=0 


(3.51) 


(3.52) 


At large n one can shift n 
1 


Ip T(p)(-b 0 )P 


n — p to obtain 
E Moa(Q 2 ))> + l)! 


Large n 


T(n + 3 — p) 


_n(n + l)r(n + 1 — p)_ 


(3.53) 


One can see that the expression between square brackets goes to unity at large n, 
independent of p. The factor 1/T{p) out front affects the residue but not the position of 
the singularity on the Borel plane at — l/&o- hr fact one does not expect any suppression 
from large values of p in the final contribution, as the coefficient accompanying the 


Q: eff(^) P t erm a priori may behave as pi [20, 211. This effectively renders the residue of 
the UV renormalon incalculable. 

Our discussion has been simplified by not including the effect of anomalous dimensions of 
the operators Oi of eq. ( |3.49| ). They introduce contributions such as [a(Q 2 ) log {k 2 /Q 2 )\ m 
for an arbitrary power m in the integral ( p.50| ). Taking into account anomalous dimen¬ 
sions makes the analysis more cumbersome but does not change the main conclusion 
about the dominance of four-fermion operators (see refs. R,|2l| for details). 


What message does one learn from these calculations? 

For the class of QCD contributions which are obtained with the “amputated” generating 


functional of eq. (T8) we argue that, when the momentum transfer in the effective 
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charges exchanged are very large, only two types of effective local operators can appear, 
depending on whether the large momentum flows to the external vertex or not. The 
local operators are respectively vertex-like operators which connect two fermions to the 
relevant external source (like the case illustrated in Fig. 4,) and four-fermion box-like 
operators (like the case illustrated in Fig. 5.) As we have seen the contribution from 

the dimension d = 6 four-fermion operators dominates because of the two powers of 
k 2 

log factors from each fermion loop. In the large N c limit, the only d = 6 four-fermion 
operators which are allowed by the chiral symmetry properties of perturbation theory 
are the two operators of the ENJL model. Four-fermion tensor operators like 


E 


ij=flavour 




'if 


R 


’ R a W 


(3.54) 


which is chirally invariant and where colour is summed over within the parentheses, 
can be seen to vanish identically for instance by employing the identity a^ u ( 1 ± 75 ) = 

In theories like QED, where the UV renornralon can be interpreted as the signal that 
non-perturbative contributions must exist, there are indications from lattice simulations 
that four-fermion operators may indeed play an important role in non perturbative 
dynamics |2(|. In QCD it is difficult to imagine how any non-perturbative dynamics may 
emerge from the study of UV renormalons since they involve only very high momenta. 
Our point of view is that the richness revealed by the UV renormalon structure may 
perhaps be used as well in QCD to learn about generic features of non-perturbative 
physics provided one studies a regime where there is some infrared momenta involved. 
The only possible choice in the two-point functions which we are studying seems to be 
the regime where the external momenta Q 2 is taken to be very small instead of very 
large as is usually assumed. One can think of many physical processes which involve 
integrals over the external field Euclidean momenta Q 2 of these two-point functions all 
the way down to zero. The hadronic vacuum polarization contribution to the anomalous 
magnetic moment of the muon is a well known example. (See e.g. ref. 47] and references 
therein.) It is this “atypical” situation of small external momenta which we next want 
to explore. 


9 We thank J. Bijnens for pointing this out to us. 
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4 The fate of UV—Renormalons in Two—Point Functions at 
low Q 2 values. 

We wish now to discuss what happens to the UV renormalon contributions to two-point 
functions when the external Euclidean momenta Q 2 becomes smaller and smaller. We shall 
take as a starting point the contribution of two chains to the vector two-point function in 
eq. (3.15). One can readily obtain from it the following leading contribution to the Adler 
function 


ntf)= 


167T 2 


N c 1 \ 




2 2-7T / 


1 


9 


f°° dk | 
IQ 2 


“eff ^e) 


log 2 {k 2 E /Q 2 ) 


(4.1) 


This integral is well defined as long as Q 2 > A 2 ; however, for Q 2 < A 2 there is a pole 
in the integration range. Let us examine this more closely. With the change of variables 
w = a(fi)bo\og(k 2 E /Q 2 ), we can rewrite eq. ( 3.16 ) as follows 


A{Q 2 ) = 


N r . 


N r . 1 


167T 2 V 2 27t6o/ \9/ 6 oa(^ 2 ) jo 


1 


dw e 6 o q (m 2 ) 


w 


(1 + wq + w) 2 


(4.2) 


where a(// 2 ) is the running coupling constant defined at a scale // 2 > A 2 , where a perturbative 
expansion in powers of ct(/x 2 ) makes sense a priori. In this equation 


wq = a(n)b 0 log {Q 2 /^ 2 ), 


(4.3) 


and in the range Q 2 < h? L one finds that wq < — 1 so that, indeed, the above integrand has a 
double pole at w = — 1 — wq, i.e. within the integration range J^] Consequently this integral 
will be ambiguous depending on how one goes around the double pole. This ambiguity turns 
out to be //-independent and reads 


M«?V<»! = - (f 


27T&0 


2 Q 2 (l 2 A 2 A 2N 

~PT E lo r - 21 °g7w 


A| 


9 


Q- 


Qr 


(4.4) 


where 1C is an unknown constant parameterizing the ambiguity and of 0(N®). 

We think this is quite a remarkable result! The ambiguity turns out to be of the same type 
as the insertion in the vector current two-point function n(Q 2 ) of the local d = 6 composite 
operator 

87 t 2 Gv 


o = — 


a,b=flavour 


E (^Ql) + (L - R) 


(4.5) 


provided one interprets h? L as the momentum cutoff in the loops and A x ~ A^. This is one 
of the four-fermion operators appearing in the ENJL model [28,29| with the identification 


4GV = -1C 


1 N r 


v 27rfeo 2 

The floating scale k E in eq. (3-13) has now turned into a hard scale A 2 . 


(4.6) 


J One can easily check that if Q 2 > Af one may choose fi = Q 2 and integrate by parts to recover eq. (3.16). 
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We can repeat the same calculation starting from the integral in eq. (O for the contri¬ 
bution of the one-chain renormalon. The ambiguity is again /U-independent but now it has 

A 2 

one power less of the log jjf. 


6U(Q 2 ) 


i one chain 

Iq 2 <a 2 ~ 


N 2 Q 2 , A 2 


b 0 A 2 


log 


Q 2 


(4.7) 


and consequently this ambiguity is subleading with respect to that of eq. ([O]) in the low-Q 2 
limit. In Section [3] we have seen that, for Q 2 large, the two-chain renormalon dominates over 
the one-chain one at large orders of perturbation theory . What we see now is that, for Q 2 
small, a similar dominance effect appears as well, and it manifests itself in the form of an 
extra log in the order of the ambiguity . The ambiguity in (4.7) can also be reproduced 
by the insertion of a local d = 6 operator, which in this case is 


N c 
bo A| 


q(x)rd a F a/1 (x)q(x ) 


(4.8) 


with F afl (x) being the field strength tensor of the external vector field v fl (x). The normaliza¬ 
tion out front is of 0(N®) but is ambiguous. The operator in eq. ([4.8|) again coincides with 


there is a constant scale in the denominator 


the effective operator at the scale Jze of eq. (343) but with the important difference that now 
re is a const an 
If one defines 


a{Q 2 ) = 


) 


(l + boa(fi 2 ) log(Q 2 //r 2 )) 


(4.9) 


for Q 2 < A 2 l , i.e. if one uses analytic continuation, expressions such as eq. (|4.2|) acquire a 
more familiar form. Notice that the a(Q 2 ) defined in this way becomes negative for Q 2 < A 2 L . 
By making the change of variables 


/ ol{Q~ 

w = — w 


a(/r 2 


one can check that eq. (|4.2| ) becomes 

A(Q 2 ) = - 


N 3 

1 V C 


9(47t) 4 6q b 0 | a(Q 2 


poo _ w' njJ 2 

/ dw 1 e 6 ol“W 2 )l --—r- , 

/ o (l-w') z 


(4.10) 


(4.11) 


and the relevant expansion parameter is now | a(Q 2 ) |, the absolute value of a(Q 2 ). In the 
ordinary situation where Q 2 > A 2 one can follow the same steps and show that a plus sign 
appears instead in the denominator, i.e. it is of the form (1 + u/) 2 , as it should be. Of course 
in this case eq. ([O]) is not merely a definition but the true solution of the renormalization 
group equation. 


We shall next argue that it is not so surprising to find that UV renormalons become real 
singularities when Q 2 < A 2 . It is helpful for this purpose to look at the Q 2 ~k^j plot in Fig. 8. 
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Fig. 8 Integration regions of the virtual Euclidean momentum K E for a fixed Euclidean 
external momenta Q 2 . 


The integration regions 0 < k E < Q 2 and Q 2 < k E < oo shown in Fig. 8 are at the origin of 
the appearance of IR renormalons and UV renormalons in two-point functions, respectively. 
(Fig. 8 corresponds to the precise case of the one-chain renormalon, but it can be easily 
generalized to an arbitrary number of chains by introducing one more dimension for each 
renormalon chain.) For reference, we also show in this plot the lines Q 2 = and k 2 E = A|. 
In the conventional study of IR renormalons, the external Euclidean momenta is always chosen 
to be Q 2 A|; however, regardless of how large Q 2 is taken, there will always be a region in 
the virtual k E integration which is below and which leads to IR renormalon poles in the 
Borel plane. As discussed in Section || their appearance is welcome because they reflect the 
limitations of perturbation theory and indicate the presence of non-perturbative -power 
corrections, as indeed the OPE in the physical vacuum suggests. What we are proposing here 
is another way to approach the non perturbative regime of QCD. We want to explore the 
kind of ambiguities which appear when the external Q 2 is chosen to be below h? L . Then, with 
Q 2 < A|, there is always a region in the virtual Q 2 < k E < oo integration region which is 
also below A|, (the darker triangle in Fig. 8.) It is precisely this integration region which 
is responsible for the promotion of the UV renormalons to poles on the right hand side of 
the Borel plane. From considering IR renormalons we are learning what type of corrections 
to the perturbative evaluation of two-point functions to expect, coming from the large Q 2 
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asymptotic regime. In the case of UV renormalons we are learning what type of effective local 
operators govern the low-energy physics of the same two-point function, when the virtual 
Euclidean momenta kg is integrated all the way from oo down to the low-Q 2 regime . The 
perturbative approach to these two non-perturbative limits is very different, and it is therefore 
not surprising that the physics in the two cases is also different. We insist on the fact that 
both limits involve an integration region of virtual momenta below the Landau scale; and in 
that respect there is no reason to expect one limit to be more physical than the other. 

From another point of view let us notice that eq. ( |4.9| ) implies that a(Q 2 ) —> 0~ as Q 2 —> 0. 
This means that the behaviour of a(Q 2 ) extracted by analytic continuation in the low-energy 
region is governed by the trivial IR fixed point at a(Q 2 ) = 0~. In this sense it is pretty much 
like an abelian theory, only that the expansion parameter is | a(Q 2 ) |= — a(Q 2 ) > 0 ; (see 
eq. ( fhTTD ) 0- We know that an abelian theory is not asymptotically free and consequently 
that its UV renormalons sit on the right hand side of the Borel plane. This is what happens 
here as well. It is then reasonable that, as seen from the extreme low-energy end, there may 
be local higher dimensional operators suppressed by the large scale since it is the only 
scale available. As the energy goes up and eventually goes over this scale A l, the dynamics 
is ultimately governed by the true UV fixed point of QCD at a(Q 2 ) = 0 + , as follows from 
its nonabelian character. As far as renormalons are concerned, if we had to exemplify this, a 
similar situation could be like having QED embedded in a nonabelian GUT model like, e.g., 
SU( 5). 


We turn next to the discussion of the pseudoscalar two-point function in the same low-Q 2 
regime. The analysis runs parallel to the vector one. As far as the two-chain renormalon is 
concerned one sees that the leading result eq. (|3.43| ) is proportional to the one in the vector 
channel, eq. (3.15). One therefore concludes that the ambiguity is 


5P(Q 2 


itwo chains __ ir A c 

'Q 2<k % ~ ^ 16vr 2 


— ) 2 
2 irboJ 


91 

A! 





(4.12) 


with /C the same constant appearing in eq. 0- This ambiguity can also be reproduced 
by a local four-fermion operator; an operator exactly of the form of the scalar-pseudoscalar 
four-fermion operator which appears in the ENJL model: 


° = 8 1 (ql( x )qU x )) ■ 

c x a,b=flavour 


The identification with the resulting perturbative ambiguity requires now 


G s = 


4/C 


f 1 N c 
\27r6o 2 


2 


(4.13) 


(4.14) 


Assuming that the ambiguity 1C is a measure of the final non-perturbative value for 
the coefficients of the ENJL-like four-fermion operators would lead to the conclusion that 
Gs/Gy = 16. However, one should not forget that this result has been obtained at the 
two-chain renormalon level. As we shall next show, higher number of chains also contribute 

n The one-loop coupling boa(Q 2 ) = (logQ 2 /Aj,) 1 is invariant under the change a(Q 2 ) —> —a(Q 2 ) , 
b 0 —> — b 0 . 
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with the same Q 2 /dependence as in eqs. (4.4,4.12), with the result that not even the ratio 
Gs/Gy can be fixed from the lowest non-trivial two-chain calculation . 


When we examined UV renormalons in the high-Q 2 regime in 3.2.4 we saw that diagrams 


with an increasing number of chains are not suppressed with respect to the contribution 
coming from two chains. This feature persists when Q 2 < A 2 . Let us consider the integral in 
(3.5C) again, but now for Q 2 < A|: 

4 

Q 2 ’ 


I P = N C Q 2 


IQ 2 


dk 2 E 


N cU e t f(fc|) log : 


(4.15) 


where we have introduced convenient factors of N c to match the dependence of a leading 
contribution at large N c . Defining w = — a(Q 2 )bo log(k E /Q 2 ), where a(Q 2 ) < 0 is given by 


equation (L9), one finds 


I P = N C 


[N c a(Q 2 )Y 


dw e (- b o a (Q~)) 
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(1 - w) p 


[-b 0 a(Q 2 )] 3 

This integral has a pole at w = 1 of multiplicity p. Therefore the ambiguity is given by 

[nmq 2 )] p 


(4.16) 


5 I p = /C N c 
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p (~ b o a (Q 2 )) qn 2 

\dw) 



, o - - ■ - . (4.17) 

[—6 0 a(Q 2 )] 3 ' J - ' 

Each time the derivative d/dw acts on the exponential it brings down a factor (&ock(<3 2 )) -1 , 
hence a power of \og(Q 2 /A? L ). When Q 2 —► 0, the leading contribution is given by 


6I P = 


-1 


K Nr, 


where 


(p- 1)! 

stands for subleading terms when Q 2 


—\ P 9l (log 2 ^ + 

bo) A 2 L { S Q 2 + 


(4.18) 


<< 


. Consequently one obtains once 
more the Q 2 dependence of eqs. (4.4) and (4. 3 ) □_ and therefore the ambiguity can be 
reproduced by the insertion of the ENJL operators ( |4.5| ) and ( 4. 13| ) , although with different 
coefficients Gs , Gy. The calculation of these coefficients for an infinite number of chains looks 
like a much harder problem since there is no obvious parameter that dictates that a higher 
number of chains should give a smaller contribution. Therefore there is no reason to believe 
the relation Gs = 16GV will be preserved beyond two chains and, for instance, the Gs and 
Gy corresponding to the three-chain double-box diagram come in the ratio Gs/Gy = 64. 
From a practical point of view, one may just as well leave the coefficients Gs and Gy as free 


phenomenological parameters to be fitted to some experimental input [28, 291. This is just 
the usual situation. Renormalons help one to make good guesses for parameterizing non- 
perturbative physics, but they are not capable of producing quantitative results: in our case, 
Gs and Gy appear to be incalculable quantities. 

There is another point of view, concerning the ratio Gs/Gy, which has been advocated by 
the authors of ref. [ 481. Their claim is that the consistency of the successes of the QCD sum 
rules with the ENJL model requires Gs Gy. If that was the case, one may then perhaps 
be tempted to take more seriously the quantitative perturbation theory results above. The 
analysis in ref. |2f|, however, shows that the success of the ENJL model in predicting the 
0(p 4 ) Li coupling constants of the chiral Lagrangian hinges on the fact that Gy has to be as 
large as Gs 0- A detailed phenomenological analysis of this issue is under investigation [49]. 


12 This persists even after use of a two-loop (3 function. See Appendix A. 

13 It is in particular the constants L 5 and Lg what requires the axial coupling constant of the constituent 


quarks <ja to be — 0.6 and hence Gy ~ Gs- 


29 
































5 Conclusions and Outlook. 


The analysis of ultraviolet renormalons, in the way we just described, seems helpful in bridging 
the gap between high energy and low energy in QCD. It offers an arena where it begins to be 
possible to discuss generic dynamical patterns of the effective low energy chiral Lagrangian of 
QCD. At high energies UV renormalons hint at the emergence of dimension d = 6 four-fermion 
operators from the large orders of perturbation theory [2^, 2l|]. What we have seen here is 
that, at low energies, they may even “feel” the existence of the scale A x in the normalization 
of these four-fermion operators. There is a correspondence between the leading n! behaviour 
at large orders of perturbation theory and the leading Q 2 —* 0 behaviour as both are due to 
the insertion of the same effective operator. 

Perturbation theory (to all orders) and analyticity in the coupling constant has led us 
this far. Let us now turn to some physical comments. The picture which emerges is obvi¬ 
ously incomplete as, so far, there is no signal of spontaneous chiral symmetry breaking or 
confinement. Neither of these features should be expected however from perturbation theory 
arguments alone, which after all is the framework of renormalon calculus. As we have seen, 
the analysis of UV renormalons in the regime of low external Q 2 momenta provides a plausible 
link with the ENJL model. This model however is only supposed to be meaningful in the 
event of a clear separation between the confinement scale A con f and the spontaneous chiral 
symmetry scale A x , as it is an effective Lagrangian description in between the two scales with 
A con f << A x . That these two scales may be separated widely enough seems to be backed 
up by the phenomenological success of the constituent chiral quark model |33|] , which follows 
naturally from the ENJL model, and wherein there is a succesful description in terms of “con¬ 
stituent massive quarks”. For this to be possible within the UV renormalon approach which 
has been advocated above, the coupling constant Gs in the d = 6 four-fermion operator ( 4.13 ) 
has to be larger than unity. It is very difficult to justify why Gs should have this size within 
perturbation theory alone. However we know that there is a (rather similar!) case —the gluon 
condensate— where the ambiguities foreseen in perturbation theory via the IR renormalon 
analysis are finally realized by a rather large vacuum expectation value of non-perturbative 
origin. There also, renormalons only signal the appearance of non-perturbative contributions 
in the form of Q~ 4 terms, and it is then a matter of non-perturbative dynamics what finally 
makes the numerator large. It is not unreasonable to assume that the coupling constants Gs 
and Gy of the four-fermion operators which govern the leading UV renormalon effects may 
behave similarly to the gluon condensate and turn out to have large values as well in the real 
world. 

We have seen that, at least within the framework of the “amputated” effective action in 
(|2 . 8| ) , four-fermion operators govern the properties of the leading UV renormalon. Assuming 
that, eventually beyond perturbation theory, chiral symmetry is spontaneously broken by the 
dimension-six scalar four fermion operator with Gs > 1, the question one may ask is: can 
one safely neglect operators with dimensions higher than six? Although they are associated 
with subleading renormalons in perturbation theory; beyond perturbation theory the counting 
of dimensions changes for the dimension-six scalar four-fermion operator with Gs > 1. Its 
effects become suppressed only by a 1/log(A 2 /Q 2 ) 2 instead of Q 2 / A 2 . This enhancement 
partially takes place also in scalar four-fermion operators with higher dimensions. The net 
effect there, however, is that these only modify somewhat the relationship between the mass 
of the scalar particle and the mass of the constituent quark (Ms ^ 2 Mq) but they do 
not change the spectrum or the interactions [50]. In the vector channel, higher-than-six 
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dimensional operators only modify the results of the ENJL model at values of Q 2 Ml 
but not at lower energy. 

More work than what we have presented here is clearly needed and seems worthwhile. 
What is at stake is that it may be finally possible to understand why the ENJL model, 
and hence the constituent chiral quark model, does so well in describing low-energy QCD 
phenomenology [^, 291. We think that the patterns unraveled so far must not be just a 
coincidence. 
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Appendix A 

In this appendix we shall outline how to generalize our previous results to the case of a 
two-coefficient (3 function such as 


a(p z )/3(a(p z )) = —boa(p 2 j 2 — ha(ii z ) 6 . 


,2\3 


(A.l) 


The use of this (3 function is interesting since i) it has a finite radius of convergence in a [ 531 (it 
is a polynomial), and ii) any /3 function can be brought to this form by making a perturbative 
redefinition of the coupling constant jT|. 

The equation 
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Notice that equation ( |A.3| ) implies that 
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The integral we want to study is 
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for n a negative integer with n < —1 and p a natural number, p > 2. 


Making the change of variables [ 521 
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where z n = n/bo < 0 and Q 2 « A|, so a(Q 2 ) ~ (bo log(Q 2 /A 2 )) 1 , one obtains 


I = 


dz e 
—n 


/ \ 1+^n 

h-i) 


1 


[a(Q 2 


z n 

z 


+ h} 11 ~ p 


oc(Q 2 ) b 0 J_ 

zb\ h / z 
- -To - 72 log 1 - — 


. z n b 0 a(Q 2 ) z n b 2 o b% 


z 

Zn / J 


(A.6) 


(A.7) 


(A.8) 
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where 5 n = nb\/b q < 0. 

Since z n < 0, the integration region contains two singularities that, in the limit a(Q 2 ) 
CD, collapse to a branch point at z = z n . Therefore 


SI = aiQ 2 )?- 1 [ 

Jo 


°° dz 
o ~n 


Disc 


FIT" 


zb\ b\ 


hence 


51 ~ a(Q 2 ) p ~ 1 I —Disc 


dz 


/ o —n 

which can still be approximated as 


z n b 0 a(Q 2 ) z n b 2 0 
e ~ Z (^(Q T ) + b5 

b-t) 


- 79 lo S 1 - 


n 2 


(A.9) 


P+Sn 


z n b 0 a(Q 2 )_ 


(A.10) 


si ~ r~ 2i e -<T?y + fe) Disc / b _ D) 

&0 Jo -n IV -W 


-P-Sn 


(A.ll) 


Given that Disc ^1 — j-'j 1 = JC( 1 — e * 27n5n ) | 1 — z/z n p <5n . where 1C is an arbitrary 

constant, one obtains 


£( 1 ~ e ") e —Znbi/bo 

—n(bi/bo) 5n n 


r( 1 -p-« (A) log 2 ^ + 


(A-12) 


where the ellipses stand for terms subleading in the limit Q 2 /h? L —> 0. One can check that 
the —> 0 limit of the previous expression is equivalent to eq. ( |4.18 ) in the text, which was 
obtained in the case b\ = 0. 

One finally sees from eq. ( |A.12| ) that the leading contribution (n = —1) produces again 
an ambiguity ~ (Q 2 /A 2 ) log 2 (A| / /Q 2 ), exactly as in the b\ = 0 case discussed in the text. 
Therefore our conclusions are also valid in the general case of the (3 function (AT). 


Appendix B 


We would like to comment here on the issue of the “freezing” of the strong coupling constant 
at low energy. 

Sometimes one finds in the literature [31] proposals for a non-perturbative behaviour of 
the coupling constant in the low-energy region of the form 


a(Q 2 ) = 


1 


bo log Q2 ^? 2 


(B.l) 


This form clearly guarantees the right high-Q 2 behaviour dictated by perturbation theory. It 
can be obtained from the following RG equation: 


Q 2 


da(Q 2 ) 
dQ 2 


= a{Q 2 )(3(a{Q 2 )) . 


(B.2) 
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with 


/3(«(Q 2 )) 


b 0 d(Q 2 ) 


C 2 i 
1 - e *> 0 &(q‘ 2 ) 

A 2 


(B.3) 


The constant C is supposed to encode the non-perturbative dynamics. As C 2 —> 0 one 
recovers the usual situation in perturbation theory at one loop. In this equation C 2 > A 2 > 0 
so that a(Q 2 ) does not have a Landau pole. 

What we would like to point out here is that chiral symmetry breaking is a strong con¬ 
straint for this proposal, at least in the simplest form of eq. ( [B.1D , as we now explain. 

If we take the Adler function A(Q 2 ) (defined in the main text) one knows from general 
arguments of spontaneous chiral symmetry breaking in the large- N c limit that A(Q 2 ) ~ Q 2 
as Q 2 —> 0 Q and in particular .4.(0) = 0. 

In perturbation theory (i.e. at large Q 2 ) the Adler function is given by 


A (Q 2 ) 


1 

87r 2 



(B.4) 


where a(Q 2 ) is the usual perturbative coupling constant at one loop, say. If a(Q 2 ) is supposed 
to incorporate all the non-perturbative dynamics by deviating from a(Q 2 ) as Q 2 —> 0 one 
would expect that 


A (Q 2 ) 


1 

87T 2 



(B.5) 


would be the natural answer. However we saw that freezing leads to a(Q 2 ) 
s Q 2 —> 0 and therefore .4(0) A 0, in conflict with chiral symmetry. 


constant (> 0) 


The problem seems rather generic for coupling constants obeying an equation such as (B.2) 
for if .4(0) has to vanish, d(Q 2 ) has to become negative at low Q 2 . But since d (Q 2 ) is positive 
at high Q 2 it follows, if a(Q 2 ) is a continuous function of its variable, that there must be an 
intermediate point at which d(Q 2 ) is small and its slope is positive, which is impossible since 
its slope is the /3 function which, for small a(Q 2 ), has to be negative according to asymptotic 
freedom. 

We would like to emphasize that the argument we presented above is not intended as a 
“NO-GO” theorem but as a sort of illustrative warning signal. One way out is to invalidate 
eq. ( |B.5| ), because more powers of a(Q 2 ) cannot be neglected or for any other reason. For 
instance, there are models that are successful in describing diffractive phenomena in which it 
argued that the freezing coupling constant QB.1|) also comes with a Q 2 -dependent dynamical 


is 


mass for the gluon. In this case eq. ( |B.5| ) does not hold and one may then be safe. As it 
turns out, however, in this case it takes a more detailed analysis of the Schwinger-Dyson 
equations to reveal that, actually, these models do not seem capable of generating enough 
chiral symmetry breaking in the end [f)4| . 

Another possibility is that the freezing of the coupling constant may have a more sophisti¬ 
cated dynamical origin. For instance, if one allows that the RG equation be governed by a (3 
function that depends on both d(Q 2 ) and Q 2 , /3^d(Q 2 ), Q 2 ^j , then eq. (BA) may as a matter 


14 In the large- N c limit only one-meson states contribute to A(Q 2 ). Since these states are not massless in 
the chiral limit, dimensional analysis imposes that A(Q 2 ) ~ Q 2 /M 2 with M the corresponding meson mass. 
15 As a matter of fact d(0)/7r < 0.3 |3lj; i.e. a relatively small correction to unity. 
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of fact reproduce .4(0) = 0. Just as an existing proof, let us quote one such a(Q 2 ): 


d(Q 2 ) = 


—Tie 


Q 

'IP 


A 

■q 2 


+ 


bo log (a + 


(B.6) 


for arbitrary parameters A 2 and a > 1. 

Therefore our conclusion is that spontaneous chiral symmetry breaking imposes rather 
severe restrictions on the idea of freezing. One should check that enough symmetry breaking 
can be produced before any argument relying on a particular freezing coupling constant is 
put forward. 
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